IAES International Journal of Artificial Intelligence (IJ-AI) 
Vol. 10, No. 4, December 2021, pp. 910~918 
ISSN: 2252-8938, DOI: 10.1159 1/ijai.v10.i4.pp9 10-918 0 910 


A hybrid approach to multi-depot multiple traveling salesman 
problem based on firefly algorithm and ant colony optimization 


Olief Imandira Ratu Farisi!, Budi Setiyono?, R. Imbang Danandjojo* 
‘Department of Mathematics Education, University of Nurul Jadid, Probolinggo, Indonesia 
*Department of Mathematics, Institut Teknologi Sepuluh Nopember, Surabaya, Indonesia 
Research and Development Agency of Transportation, Ministry of Transportation Republic of Indonesia, Jakarta, 


Indonesia 


Article Info 


ABSTRACT 


Article history: 


Received Feb 20, 2021 
Revised Aug 25, 2021 
Accepted Sep 12, 2021 


This study proposed a hybrid approach of firefly algorithm (FA) and ant 
colony optimization (ACO) for solving multi-depot multiple traveling 
salesman problem, a TSP with more than one salesman and departure city. 
The FA is fast converging but easily trapped into the local optimum. The 
ACO has a great ability to search for the solution but it converges slowly. To 


get a better result and convergence time, we integrate FA to find the local 


solutions and ACO to find a global solution. The local solutions of the FA are 
normalized then initialized to the quantity of pheromones for running the 
ACO. Furthermore, we experimented with the best parameters in order to 
optimize the solution. In justification, we used the sea transportation route in 
Indonesia as a case study. The experimental results showed that the hybrid 
approach of FA and ACO has superior performance with an average 
computational time of 26.90% and converges 32.75% faster than ACO. 
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1. INTRODUCTION 

Multi-depot multiple traveling salesman problem (MmTSP) is the shortest route problem which is 
the further development of traveling salesman problem (TSP). In MmTSP, more than one salesman must visit 
a set of cities exactly once and depart from different departure cities, called depots, and must return to their 
respective departure cities [1]. Finding the exact solution of the shortest route is more difficult as the number 
of visited cities increase [2]. 

Unlike TSP which has been extensively researched recently years, the MmTSP is not intensively 
studied. In fact, many applications which are related to routing problems with multi-depot can be modeled as 
MmTSP, such as the logistic distribution problem [3], vehicle routing problem [4], vehicle scheduling 
problem [5], split-delivery vehicle routing problem [6], and the location routing problem [7], [8]. Since TSP 
is NP-hard and the MmTSP includes the TSP as a particular case, so it has a higher complexity than TSP [9]. 
Therefore, it would be very useful to design an algorithm that has a good performance and produce a better 
solution of MmTSP. 

Several researchers have conducted studies on the heuristic algorithms for solving MmTSP, such as 
multiple depots and closed paths (MDCP) method [10], genetic algorithm (GA) [11], firefly algorithm (FA) 
[12], and ant colony optimization (ACO) [13], [14]. According to Ghafurian and Javadian [13], the ACO is 
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compared to Lingo 8.0 and showed that the ACO obtained an optimal solution faster than Lingo 8.0. The 
ACO is an optimization algorithm that adopts the behavior of ant colonies which was developed by Dorigo to 
solve TSP [15]. The ACO is robust and has a great ability in finding NP-hard problems solution [16]. It is 
also easily applied and understandable [17]. However, the change in the ACO’s probability distribution in 
each iteration causes the uncertainty of the ACO’s convergence [18]. Therefore, a few researchers developed 
a hybrid approach to improve the robustness of the ACO method, such as the integration of FA and ACO. 

The FA is an algorithm based on the flashing behaviour of fireflies and was proposed by Yang [19]. 
The FA is easy to implement and has a fast convergence capability but it can be trapped into local optimum 
solution [20], [21]. A combination of FA and ACO approaches has been developed in solving various 
problems, such as economic power dispatch [22], traveling salesman problem (TSP) [23], and constrained 
engineering design problems [24]. According to the researches, a hybrid approach of FA-ACO has better 
performance than ACO itself. 

In this regard to improve the accuracy and convergence time of ACO in finding the solution of 
MmTSP, we propose a hybrid approach of FA and ACO to optimize sea transportation routes in Indonesia to 
increase the effectiveness and efficiency of sea transportation services by the ministry of transportation. This 
problem can be modeled as a MmTSP since 51 ports must be visited with 5 ports being the main focus of the 
sea highway development as depots. Using the fast converge of FA, we find the local solutions with a shorter time. 
The local solutions are used for pheromone's initialization in the ACO to find the optimal solution. 


2. PROPOSED METHOD 

In this paper, we develop a hybrid FA-ACO to solve MmTSP. In this method, FA is used to find the 
local solutions then ACO will find the global solution based on FA's solutions. By finding the best 
combination of the lower possible number of ants and fireflies as well as all the parameters needed to run FA 
and ACO, we can optimize the solution and minimize the convergence time. In the ACO method, the initial 
pheromones trail is distributed equally for each edge. But, in the hybrid FA-ACO, we set the initial trail 
based on the local solutions obtained by FA. The pheromones will be added on the edge which is part of the 
best local solutions. In this way, the probability of ants visiting the edge is bigger. The stages of the proposed 
method are shown in Figure 1. 


MmTSP’s input 7’ Local Search by FA 

: 
Tahal oanonot Solution Normalization 
pheromone : | 
se Initialization of 
we pheromone based on 
Global Search by ACO \ FA’s solutions 
MmTSP’s solution 


Figure 1. Stages of hybrid firefly algorithm-ant colony optimization 


2.1. Local search by FA 

Firefly Algorithm is an algorithm that is inspired by light-based fireflies’ movements. It was first 
introduced by Yang [19] to solve continuous problems. To solve the permutation problem like MmTSP, the 
FA is discretized [25]. 

Fireflies move towards other fireflies that have a brighter light. The intensity of the light is 
influenced by the distance and absorption of light by air. The farther the distance, the smaller the intensity of 
the light. Also, the greater the absorption of light, the smaller the intensity of the light produced. 

Based on [26] work, in this study, the fireflies represent permutation solutions of MmTSP. For 
example, Figure 2 shows the permutation representation of MmTSP with city-1 and city-5 as depots. The tour 
for the first salesman is 1-8-4-10-6-1 while the second salesman's tour is 5-7-9-3-2-5. 

Since the objective function of the MmTSP is minimizing the cost, the light intensity of a firefly (J) 
is calculated as the inverse of the total cost. The smaller the total cost produced, the brighter the light intensity 
emitted. After the fireflies have their brightness, the distance between two fireflies is determined by (1), 
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Ad(i, j) 
= —— xX 


di, J) = 10 (1) 


where Ad(i, j) is the number of different edges between the permutation representation of firefly-i and 


firefly-j, and N is the number of cities. In Figure 3, three edges 4-10,10-6,7-9 in firefly-i do not belong in 
firefly-j. In other words, there are three edges 4-6, 7-10, 10-9 in firefly-j do not exist in firefly-i. Therefore, 
the distance between firefly-i and firefly-j is 3. A firefly moves toward the most attractive firefly based on the 
brightness and the distance. The attractiveness of firefly-j seen by firefly-i is any monotonic decreasing 
function in (2), 


E(i,j)=ée" (2) 


where d is the distance between the two fireflies, ¢o is the light intensity of firefly-j at zero distance, and 7 € 


[0, 2) is light absorption. If y > 0, then (i,j) = &. It means that at any distance, the attractiveness of 
firefly-j when seen by firefly-i is the same. If y — 00, then firefly-j cannot be seen at any distance because the 
value of light absorption is too large. Hence, the value of z must be tested because it will affect on the 


attractiveness of fireflies. The coefficient y is in the interval [0.01, 0.15]. 

A firefly moves toward the other firefly that has the most attractiveness. If there are no fireflies with 
greater attractiveness, it moves randomly. The number of movement's steps that must be taken by firefly-i 
which attracted to the brighter firefly-j is (3); 


x; = random(1, Ad(i, j)) (3) 


In MmTSP, the movement of firefly is the inversion mutation of its permutation representation. 
Inverse mutations are carried out in visited cities excluding depots. By using inversion mutation, the previous 
permutation represented by the firefly is not deformed. Figure 4 illustrate the movement of firefly-i. 

This movement will make changes to the solutions represented by fireflies. Since each firefly moves 
for m times, it will produce m new solutions. If there are p moving fireflies, there will be pm new solutions. 
Choose the best solutions p of p + pm existing solutions as the new population of fireflies for the next iteration. 


Pris [4 if 6 Prpst7{9 1312s) City 
i 2 3 4 . 6. o- Ss 


Order 


Figure 2. Permutation representation of MmTSP’s solution 
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Figure 3. The distance of the two fireflies is the Figure 4. Inversion mutation with 2 length 
different edges between firefly-i and firefly-j movements 


2.2. Normalization and initialization 

After FA satisfied the stopping condition, there will be p best fireflies. In several cases, there is a 
possibility that some of the fireflies have the same representation of the solution. Hence, it is necessary to 
normalize the solutions. The same solutions will be considered as one solution. Then, assume that there are q 
solutions as candidate tours after normalization. From the g candidate tours, add the pheromones on the edges 
included in the tours to initialize the pheromones of ACO. The rules for adding pheromones are based on the 
order of the best solutions which is defined as (4); 
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where (i, j) is the edge of the k-th best solution of the q existing solutions. In this way, the edge that is part of 
the first-best solution will get more pheromones than the edges on the second-best solution. Thus, the most 
passed edges of the g candidate tours will get the most pheromones addition. By setting initial pheromones 
with different quantity based on the local solutions of the FA, it will minimize the number of agents and 
iterations for running the ACO. 


2.3. Global search by ACO 

ACO is an algorithm developed by Dorigo to find the shortest path which adapts the behavior of the 
ants [15]. In finding food sources, ants produce pheromones which is useful for recognizing same-sex, other 
individuals, and groups as well as for helping the reproductive process. Moving ants will release pheromones 
on the ground and create a trail. After a certain interval, pheromones will experience evaporation. 

For finding the global solution of MmTSP, place n ants at one of the depots. For each ant, the first 
element of the tabu list is equal to its depot. Set initial trail t,;(¢) for quantity of pheromones at time ¢ on 
edge (i,j) from normalization results based on the FA’s solutions. The probability of ant-k moves from city-i 
to city-j is defined as (5): 


ol“ Inj ol? 
Pij = ) Lucugltiu ©] “niu (1? (5) 
0, otherwise 


where 77;;is the visibility namely the inverse of the total distance, and U, is the set of unvisited cities by ant-k, 
with a and f are the factors that influences the pheromones and visibility, respectively. After all of the ants 
have completed a tour, the tabu list is filled with the visited cities produced by each ant. Then, the trail 
intensity is updated using the global pheromones update rule. 

The pheromones update rule involves the evaporation on all edges and pheromones’ addition on the 
ants' passed edges. It is also affected by the distance produced by the ants. The shorter the tour, the more 
pheromones the ants leave behind. Thus, the edges with lots of pheromones will be more desirable on tours in 
the next iterations. The global pheromones update rules defined as (6): 


t, (new) =(1- p)z, (old) + Ac! (6) 


where p represents the evaporation of trail, and Atk is the addition of pheromones produced by ant-k on the 
edge (i, j) is calculated by (7), 
—- if ant — k uses (i,j) 
Ati = 4dk (7) 
0, otherwise 


where Q is a constant and d; is the total distance passed of ant-k. 

The best tour produced by the ants is saved and the tabu list is emptied to save the tour for the next 
iteration. This process will be repeated until it meets a stopping condition. At the end of the method, we will 
obtain the tour with minimum cost as the solution of MmTSP. 


3. RESULT AND DISCUSSION 
3.1. Mathematical Modeling 

Based on the mathematical model of MmTSP [27], the sea transportation route problem in Indonesia 
can be modeled as follows. Let N be a finite set of ports and C;; = [c; ar (i, jEN)be the cost matrix from port-i 
to port-j. The set N is partitioned into N = DUN’ where D be the first five-port of N as depots and N’ = 
{6,7, ...,51} be the set of ports excluding depots. There is 1; ship is placed at depot-i. 

For any ships, v; is the number of visited ports on the ship’s path from depot to port-i. Constant L 
and U are the minimum and the maximum number of ports a ship may visit, respectively; thus, 1 < vy; < 
UVi > 2. In this study, we used lower bound L = 8 and upper bound U = 12. A binary variable x; is defined 
as equal to | if the ship departure from the Ath port passes the (i, j) and x;;, is equal to 0, otherwise. If xjx4 = 
1, then 8 < v; < 11 must be satisfied; 


A hybrid approach to multi-depot multiple traveling salesman problem ... (Olief IImandira Ratu Farisi) 


9140 ISSN: 2252-8938 


min [EE TC Mg i ee, ) (8) 


keD jeN' keD ieN' jeN' 
s.t 


dx, =1.keD (9) 


et Dm HLieN' (10) 


x, + yx Xi ss =0,k Ee D,jeN' (11) 


Di Xu — etme = Oke D (12) 


jew’ jen 
v, +10), x, + > x, <1LieN' (13) 
Vv. oe + (-6)>) x, $2,iEN' (14) 

x, + Dx, <LieN' (15) 
v.-v, +12) x, +10) x, SILi# fi, pe N' (16) 


From the formulation, (8) is the objective function. As shown in (9) ensures that there is exactly 1x 
ship depart from each depot. As shown in (10) ensures that each port is visited exactly once. As shown in 
(11) and (12) are constraints that represent route continuity for ports and depots, respectively. Inequality (13) 
and (14) require upper and lower bound to the number of ports in the tour, respectively. These constraints v; 
to be equal to 1, if i is the first port in a tour. Inequality (15) prohibits tours with just one port. Inequality (16) 
is sub-tour elimination constraints (SECs) that prevent possible solutions which have closed tours of ports 
with no depots as starting and ending ports. 


3.2. Experimental results 

The experiment is implemented on a 2.2 GHz Intel(R) Core(TM) i7 computer with 4.00 GB RAM 
using MATLAB. The experiment is examined for the sea transportation route. There are 51 ports with 5 ports 
as depots. In this study, there is exactly one ship in each depot and each ship may only visit a minimum of 8 
ports as lower bound and a maximum of 12 ports as the upper bound. 

We conducted an experiment to find the best parameters for running a local search by the FA and a 
global search by the ACO. The parameters play an important role in finding a solution, as shown in Figure 5. 
The parameters used are those that can reach the optimal solution in the least possible time. The minimum 
number of iterations to run a local search by FA is 100 iterations shown in Figure 5(a). For running a global 
search by the ACO, the convergence is reached after 229 iterations shown in Figure 5(a). The result is 
obtained by running the method for 1000 iterations. If no improvement of the result after 100 successive 
iterations, then we stop the iterations. In Figure 5(b) and Figure S(c), 4 fireflies and 5 movements for each 
firefly, respectively, are enough to get the optimum solution. Figure 5(d) indicated the value of light 
absorption to get an optimal solution. We get y = 0.11. 

Table 1 shows the correlation of a and # to get an optimum solution. The value tested were: a = 
{0.5, 1,2} and B = {1,2,5}. The optimal solution is reached when a@ = 1 and 8 = 5. We used these values for 
each experiment. In Figure 5(e), the best solution is obtained when p = 0.5. The minimum number of ants 
used for global search is 20 shown in Figure 5(f). 
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Figure 5. Parameter correlation with the best solution known, (a) the number of minimum iterations needed 
by FA local search and ACO global search, (b) correlation of the number of fireflies with the best solution 
known, (c) correlation of the number of movements with the best solution known, (d) correlation of light 
absorption with the best solution known, (e) correlation of coefficient with the best solution known, 

(f) correlation of the number of ants with the best solution known 


From the experiments, the parameter settings to execute hybrid FA-ACO are 
y=0.11,0=1,B=5,p=0.5,Q=100, 4 fireflies with 5 movements, 20 ants, and 100 and 229 iterations of the FA 
and ACO, respectively. As shown in Figure 6 the comparison of performance between hybrid FA-ACO and 
ACO, Figure 6(a) presents the solution of hybrid FA-ACO and ACO after running the experiment 50 times. 
The hybrid FA-ACO reached the best solutions 31 out of 50 times by achieving an optimal solution that is 
better than ACO. It showed that the hybrid FA-ACO has better performance than ACO. 
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The performance’s comparison of hybrid FA-ACO and ACO from 50 experiments is shown in 
Table 2. Our proposed method achieved the best solution known closer to the exact solution than ACO. The 
hybrid FA-ACO has an average solution better than the ACO. The average of convergence time of the 
method is also better than ACO. It can conclude that the proposed method yields a better solution and 
computationally faster than ACO. 

Figure 6(b) shows the convergence of hybrid FA-ACO and ACO to achieve the best solution. In 
getting the best solution, which is 13828 km, the hybrid FA-ACO converged after 16.86102 seconds. 
Meanwhile, the convergence of ACO whose best solution is 13882 km km reached after 25.07218 seconds 
these results proved that our proposed method can improve the performance of ACO. It can find a better 
solution and converge faster than the ACO. 


Table 1. Correlation of and coefficient with the best solution known 


a B The best solution is known (km) Running time (s) 
0.5 1.0 21038 81.4154 
0.5 2.0 17169 80.0412 
0.5 5.0 14659 83.1577 
1.0 1.0 19088 80.0881 
1.0 2.0 16072 79.0493 
1.0 5.0 14231 81.7190 
2.0 1.0 18872 79.8144 
2.0 2.0 16532 77.6645 
2.0 5.0 15361 81.1222 
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Figure 6. The comparison of performance between hybrid FA-ACO and ACO; (a) in reaching the best 
solution known from 50 experiments, (b) convergence time in reaching the best solution 
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Table 2. Comparisons of hybrid FA-ACO and ACO from 50 experiments 


Method The best solution is known (km) Average solution (km) Average time (s) 
Hybrid FA-ACO 13828 14400.28 44.42685 
ACO 13882 14463.24 60.77814 


4. CONCLUSION 

In this paper, a hybrid approach based on FA and ACO for solving multi-depot multiple traveling 
salesman problems was presented. The proposed method has three main stages: i) local search by FA; ii) 
normalization and initialization; and iii) global search by ACO. The main purpose is to find a better solution 
and minimize the time of convergence. Parameters play an important role in running the hybrid FA-ACO. 
The best parameter settings to execute hybrid FA-ACO to solve 51 ports with 5 depots and one salesman in 
each depot are y = 0.11,a = 1,6 =5,p = 0.5,Q = 100, 4 fireflies with 5 movements, 20 ants, and 100 and 
229 iterations of the FA and ACO, respectively. The experiment results demonstrated that the hybrid FA- 
ACO provides better performance than the ACO. Our method obtained the best solution better than ACO and 
achieved an average computation time of 26.90% faster than ACO. Also, to reach the best solution, hybrid 
FA-ACO converges 32.75% faster than ACO. 
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